The occurrence and human health risk assessment of total and aflatoxin B1 in selected food commodities in Bhutan

Aflatoxins are mycotoxins that contaminate staple foods globally and pose a significant health risk. To the best of our knowledge, information on the occurrence of aflatoxins in Bhutanese diets is scarce. This study aimed to estimate the aflatoxin levels in selected foodstuffs in Bhutan and determine the health risk associated with aflatoxin exposure. Ten different types of food commodities were randomly collected from farmers’ markets, shelves of supermarkets, and wholesale and retail shops from 20 districts of the country. The samples were subjected to analysis by an enzyme-linked immunosorbent assay for both total aflatoxins (B1, B2, G1 and G2) and aflatoxin B1. Among the 315 samples included, 48.81% and 79.35% were positive for total aflatoxins and aflatoxin B1, respectively. The overall mean total aflatoxin concentration was 11.49 ± 12.83 µg/kg, and that for B1 was 17.62 ± 23.99 µg/kg. The most prevalent food commodity with the highest aflatoxin contamination was chili products. In addition, the estimated daily intake and margin of exposure to aflatoxin B1 via the consumption of chili products ranged from 0.98 to 5.34 ng kg−1 bw day−1 and from 74.90 to 408.10, indicating a risk for public health. The liver cancer risk was estimated to be 0.01 and 0.007 cancers per year per 100,000 population resulting from the consumption of chili products. The present findings revealed the presence of total aflatoxins and aflatoxin B1 in the selected samples. The margin of exposure values was exorbitant, demanding a stringent public health measure. Notably, these results suggest the need for routine monitoring of aflatoxin contamination in the region and stress rigorous safety management strategies to reduce exposure.


Occurrence of aflatoxins
A total of 315 samples consisting of 10 different categories of food commodities were analyzed for the presence of AFT (B 1 , B 2 , G 1 and G 2 ), and those samples detected with > LOD of AFT were further subjected to analysis for aflatoxin B 1 (Table 2).Of these, 48.81% of the samples were positive for AFT, and 79.35% were positive for AFB 1 .The mean levels of AFT and B 1 in all the samples were 11.49 ± 12.83 µg/kg and 17.62 ± 23.99 µg/kg, respectively.The concentrations of AFT in the 95th and 50th percentiles were 41.98 and 6.84 µg/kg, respectively, while those of AFB 1 were 85.89 and 10.21 µg/kg, respectively.
Table 3 illustrates the results of the analyzed samples based upon the type of sample and the occurrence and concentration of aflatoxins.The overall concentrations of AFT and AFB 1 ranged from less than the limit of detection to 75.84 µg/kg and 116.41 µg/kg, respectively.Among the types of food samples, chili powder had the highest mean level of total aflatoxin (19.65 µg/kg), followed by nuts (14.06 µg/kg).Similarly, chili powder was the dominant food commodity, with the highest mean concentration of aflatoxin B 1 detected (28.01 µg/kg), followed by nuts (22.96 µg/kg).Although split pulses had a higher mean, only one sample was detected with a high level of B 1 , which resulted in a higher overall mean.Notably, none of the rice samples tested had > LODs for AFT; hence, these samples were not subjected to further testing for AFB 1 .
The concentration of chili products was significantly greater than that of all other food items (p value = 0.00104).Moreover, the concentration of aflatoxin B 1 correlated well with the AFT concentration, ranging highest in chili powder, nuts and dried chili, at 27.42 ± 26.49, 23.41 ± 30.82 and 11.41 ± 4.27 µg/kg, respectively.
With respect to the mean AFT concentration in the food samples from the different districts, the highest concentration was detected in Monggar (24.19 µg/kg), followed by Wangdue Phodrang (21.7 µg/kg) and Chukkha (19.88 µg/kg).The lowest values were determined for the districts of Lhuentse, Sarpang and Trashiyantse.Similarly, aflatoxin B 1 was also found to be the most abundant toxin in the abovementioned districts, at 37.  35.31 and 34.85 µg/kg, and the least abundant were from Sarpang, Tashiyantse and Gasa.However, there was no significant difference in the concentration of aflatoxins in the food samples (p value = 0.03-0.64).Figure 1 illustrates the findings.
Bhutan is divided into four regions and there is a difference in culture, food habits, ethnicity and climatic conditions.Hence, an attempt was made to report aflatoxin contamination in different regions of the country.The supplementary table shows the different levels of total aflatoxins detected in different regions of the country (see Supplementary Table S1).The highest percentage of positive samples was detected in the western region (56.84%),followed by the central region (56.67%).However, the southern region had the highest mean aflatoxin concentration (14.76 ± 16.15 µg kg −1 ), followed by the western region (12.93 ± 15.57 µg kg −1 ).One-way ANOVA was performed to compare the means of these different regions; however, no significant differences were detected in the mean level of aflatoxins.
The results from the analyzed samples were further compared with international guidelines elsewhere (CODEX, the EU and the US FDA).In total, 43.6% of the total samples had total aflatoxin levels above the European Union maximum level of 4 µg/kg, whereas 9.1% and 4.6% of the samples had total aflatoxin levels above the CODEX and US FDA maximum limits, respectively (Table 4).To ensure maximum food safety, foodstuffs were more specifically compared with the EU regulations of minimal aflatoxin contamination.The results showed that the chili powder had the highest level of aflatoxins above the regulatory limit of 94.87%, followed by dry chili (77.7%) and peanuts (60.6%).Supplementary Table S2 shows the data.

Human health risk characterization
Human health risk estimation was conducted on the most consumed food items in the country, which were chili products (chili powder and dried chili).The results are illustrated in \* MERGEFORMAT Table 5.The exposure estimates were calculated based on the estimated daily intake (EDI) of these food items.In the present study, the daily intake of chili powder and dry chili was estimated to be 2.1 g/day and 15.92 g/day, respectively, for an adult weighing an average body weight of 60 kg, which was determined from the Bhutan Living Standard Survey 2022 20 .The EDIs for chili powder at the mean and 95th percentile exposure was 0.98 and 5.34 ng/kg −1 bw day −1 , respectively.Similarly, for the dried chili, the mean exposure was 2.96, and for 95th percentile was 5.30 ng kg −1 bw day −1 .Risk characterization was conducted based on the margin of exposure (MOE) and the liver cancer approach.MOE compares exposure levels to a toxicological threshold and is typically used to compare the health risk of various contaminants to prioritize risk management efforts 21 .The lower the MOE, the greater the risk to the consumers and a value of less than 10,000 is used to indicate health risk 22 .This is measured by the benchmark dose lower confidence limit (BMDL) for a benchmark response of 10% of 0.4 µg/kg body weight per day for the incidence of hepatocellular carcinoma (HCC) in male rats following AFB 1 exposure as per the European Food Safety Authority 23 .
The range of MOEs for BMDL 10 derived from the analysis ranged from 74.90 to 408.10.From these calculations, the MOE values were less than 10,000, indicating a public health risk due to AFB 1 exposure through the consumption of chili products.Remarkably, if the 95th percentile values for large consumers are considered, there is a significant risk from exposure to aflatoxins to the general population.The incidence of HCC was estimated by the quantitative cancer risk approach and is expressed as cancer cases per 100,000 individuals per year.The values show that, for the maximum exposure scenario, the cancer risk was estimated to be 0.01 cancer cases per 100,000 people per year for the Bhutanese population.
Figure 2 shows the plot for exposure to AFB 1 and the estimate of toxicity due to consumption of dried chili and chili powder in the Bhutanese population.The intersection of the estimates of toxicity and the estimates of exposure revealed the risk of AFB 1 .A higher risk and toxicity of AFB 1 are illustrated by the closer point of intersection to the x-axis and farther from the y-axis.According to the risk 21 matrix, chili powder has a greater risk of daily exposure to AFB 1 than does dry chili.

Discussion
This study reports the first findings on the status of aflatoxin contamination in major staple foods consumed in the region and its associated human health risk.The report revealed the presence of total aflatoxins and aflatoxin B 1 in almost all of the samples included in the test.The MOE values demand a stringent public health measure.
An initial step for determining the true burden of aflatoxin contamination is the use of proper test assays 24 .For a specific test to be suitable for its intended use, the quality, reliability and consistency of the analytical results are integral.For this purpose, recovery analysis was performed to determine the actual concentration of an analyte in the samples, and we reported recoveries ranging from 87.19 to 112.02% for both AFT and AFB 1 .These percentages were within the acceptable range of the AOAC, CODEX and EU regulations of 60-120% 25 ; therefore, the analytical method was reliable and in accordance with international guidelines.
The present study revealed a positive frequency of 48.85% for AFT and 79.35% for AFB 1 .This indicates that the Bhutanese population is exposed to the carcinogenic toxins that are present in the daily diet.Similar studies performed in Thailand reported an occurrence of 38.9% 26 ; 27.1%, in Zhejiang Province, China 27 ; 24.3%, in Vietnam 28 ; 78.0%, in Nepal 29 ; and 61.3%, in India 30 .The presence of these aflatoxins is inevitable, as factors such as changes in temperature and water availability related to actual climate changes (increased temperature, heavy www.nature.com/scientificreports/rainfalls and droughts) all modulate mold growth and the production of mycotoxins 31 .This demands a safe and effective multifaceted approach for combating food and feed contamination.Although most developed countries have stringent rules for regularly monitoring these toxins in food, ensuring food insufficiency and security is a major challenge for developed and developing countries.
Studies across the region have reported AFB 1 concentrations ranging from 48-383 µg/kg in corn samples in India, 0.04-21.30µg/kg in rice in Pakistan and 1-34.80 µg/kg in corn in Vietnam 32 .In the present study, the overall AFT and AFB 1 concentrations ranged from < LOD to 75.84 µg/kg to 116.4 µg/kg for all the food commodities.This is very alarming since getting exposed to such a level of toxins leads to severe health consequences.A study by Angele et al. revealed that a concentration of 0.099 µg/kg in eggs and 0.525 µg/L in milk are risk factors that could increase the incidence and prevalence of malnutrition and cancer in Cameroon 33 .A similar study in India by Asim et al. revealed a threefold increase in the risk of hepatocellular cancer with respect to urinary AFB 1 and a higher incidence of aflatoxin contamination in the region 34 .Currently, 27.5% of the children in the country are stunted, almost half (42.6%) of the stunted children are severely shunted 35 , and the incidence of liver cancer is 5.3% 36 .The association between exposure to aflatoxins and the probability of stunting and the development of liver cancer should be considered.The question arises whether these health implications are attributed to issues related to the contamination of aflatoxin in the food that the population consumes daily.An exposure assessment study using appropriate biomarkers for aflatoxin exposure involving the recruitment of liver cancer patients and shunted children is necessary.
Notably, there was a clear difference in the concentrations of both AFT and AFB 1 among the different districts of the country.This may be attributed to the geographical locations of the districts.For instance, the Monggar district, which falls in the eastern region of the country, had an overall mean concentration of 24.19 µg/kg, and the Chukkha district, which lies in the southern region, had 19.88 µg/kg.The climate conditions vary among these regions; thus, differences were observed.A study by Battilani et al. reported that in Europe, aflatoxin B 1 contamination in maize has increased due to climate change 37 .This finding is further supported by the findings of Benigni et al. and Sarah et al., whose findings correlated well with the current findings 38,39 .
The most prevalent foods contaminated with aflatoxins reported in this study are chili products, followed by nuts.Contamination of chili powders by aflatoxins has been described by several surveys from studies in the region.In India, almost 66.7% of the samples had values above 10 µg/kg 40 ; 86.7%, in Thailand 41 ; and 100% of the samples > 20 µg/kg, in Bangladesh 42 .Similarly, in the present study, 84.6% of the chili powder samples had a concentration greater than 10 µg/kg, which is in concordance with the findings of other studies.Chili powder is produced locally using traditional methods by sun drying the chilies and grinding it to powder.In this process, there is a maximum likelihood of mold formation.The occurrence of aflatoxin contamination is determined by factors such as pre-harvest, post-harvest, storage conditions, temperature and humidity 15,32 .Therefore, prevention measures such as treating fungal contamination in crops, appropriately handling food items, ensuring effective food regulation systems and providing mass advice are essential for minimizing aflatoxin contamination in food.
The estimated daily intake of AFB 1 via the consumption of chili products was determined to be 0.98 and 2.96 ng/kg −1 bw day −1 in the present study.Studies elsewhere in Ethiopia reported an EDI of 0.13-15.8ng/kg −1 www.nature.com/scientificreports/bw day −143 , and those in Malaysia were 0.21-1.32ng/kg −1 bw day −144 .This discrepancy might be due to differences in the concentrations and consumption data collected during the particular study period.Nevertheless, the findings reveal that the risk of aflatoxin B 1 exposure originates from both the high occurrence of aflatoxin B 1 contamination and the high consumption of chili products by the population.Taken together, these findings indicate that the current regulatory limits of 20 µg/kg may still not be low enough to eliminate these concerns.Based on the current findings, it is imperative that the intake of contaminated chili products is of great importance for public health since all the MOE values are below the safe threshold value of 10,000.The MOE values of AFB 1 recorded in the present study ranged from 74.90-408.10,which are comparable to those reported for chili products from Sri Lanka ranging from 45-78 and 997-909 in Turkey 45,46 , exposing the respective consumers.

Conclusion
This study illustrates the presence of total and aflatoxin B 1 in selected food commodities from the country, providing the first insights into their occurrence and human health risk assessment.The incidence of both total and aflatoxin B 1 contamination in food samples was greater than that in half of the samples.The margin of exposure was less than the safety limit of less than 10,000, indicating a serious public health concern.In addition, the liver cancer risk was estimated to be 0.01 and 0.007 cancers per year per 100,000 population resulting from the consumption of chili products.The data presented here provide evidence that addressing aflatoxin contamination and exposure to aflatoxins is necessary and may cause severe health consequences for the population.

Study area and site
The Bhutan region is a subtropical region and has twenty administrative districts.The study site was located in each district, and samples were collected from each site.Figure 3 represents the districts of the country and the number of samples collected from each district.

Sampling
A total of 10 different categories of food commodities were randomly sampled and purchased between September and October 2022 from grocery and supermarket shelves.Each sample was homogenized by grinding, and an aliquot of 20 g was kept in zip lock bags and stored at − 20 °C until analysis.Based on the literature on the highest likelihood of aflatoxin contamination, the following types of samples were collected and subjected to chemical analysis (Table 6).

Chemicals
Commercially available detection kits (AgraQuant, Romer labs, Inc., Newark, DE, USA) for both total aflatoxin (AFB 1 , AFB 2 , AFG 1 and AFG 2 ) and AFB 1 were used for detection and quantification in the present study.The methanol used for extraction was of HPLC grade, and all the other chemicals used were of ACS grade or higher.

Sample analysis
Briefly, 20 g of representative sample was meshed and extracted with 100 ml of 70% methanol.The solution was vortexed and settled, and the top aqueous layers were filtered using a Whatman No. 1 filter.The extracted

Table 2 .
Levels of total aflatoxins and aflatoxin B 1 .
Figure 1.Burden of aflatoxins in different districts of Bhutan.

Table 4 .
Proportion of samples above the maximum regulatory limits.